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При помощи расчетов методом теории функционала плотно- 
сти систематически исследован процесс образования ком- 
плексов азотсодержащих молекул с малыми кластерами меди 
Си; (п=1-7, 13). Показано, что молекулы К М=У (У=СКВ.В., 
МЕБ›, О) перспективны для поиска присадок для кластеров 
меди, так как являются синтетически доступными, имеют 
возможность проявлять восстановительные свойства, прочно 
связываются с атомами меди и не искажают исходную струк- 
туру кластера. Используя любой объемный заместитель К, 
можно блокировать доступ к большой площади поверхности 
кластера для агрессивных соединений. Комплексы кислорода 
на поверхности медного кластера существенно уступают в 
прочности структурам, образованным молекулами Е М=У 
(У=СК.Бз, МВЬ, О). В зависимости от размера кластера сила 
взаимодействия изменяется синусоидально от минимума к 
максимуму. 
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Введение. В последнее время большой интерес вызывают исследования структуры и новые области приме- 


нения кластеров металлов группы меди [1-3]. Нанокластеры этих металлов перспективны как для использования в 
медицине в качестве маркеров или части молекулярной структуры медицинских препаратов, так и в химической тех- 
нологии (селективный катализ), оптике, электронике и ряде других областей техники и технологии [4—7]. Примером 
такого использования является применение наночастиц меди в триботехнике (раздел трибологии, науки о трении и 
изнашивании). В частности, большой интерес представляет собой использование кластеров меди в качестве присадки 
к маслам и смазкам для улучшения их триботехнических характеристик [8—10]. 

Свойства медных кластеров зависят от их размера и структуры, что в значительной мере определяется спосо- 
бами их получения и хранения. Кроме этого, важнейшей проблемой при синтезе является высокая способность нано- 
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кластеров к агрегации и укрупнению, что приводит к необходимости поиска эффективных стабилизаторов. На данный 
момент существует ряд физических, химических и биологических методов получения металлических наносистем, но 
наиболее перспективными являются методы с минимальным числом стадий и компонентов [8]. Реакции химического 
восстановления наиболее перспективны ввиду их простоты, низкой стоимости и универсальности по сравнению с 
микробиологическими и физическими методами. При этом восстановителями нередко являются органические и неор- 
ганические молекулы с активными азотными и кислородными центрами, которые и выступают в качестве стабилиза- 
торов получаемых кластерных структур [11]. Таким образом, использование азотсодержащих соединений в качестве 
стабилизатора для наночастиц меди является перспективным. 

Ранее было доказано [12], что наиболее прочными контактами между поверностью кластера меди и элемен- 
торганическими веществами отличаются соединения с азотным центром — аммиак и амины, при этом сила такого 
взаимодействия изменяется с размером кластера и при вариации среды. В связи с этим интересно изучить энергетику 
Си-М взаимодействия в зависимости от размера кластера меди и электронного строения азотного центра. 

Целью исследования было изучение влияния типичных случаев локализации электронной плотности на азот- 
ном центре на контакт с поверхностью кластеров меди в различных апротонных средах (вакуум, углеводород н- 
гептан, двуатомный спирт-этиленгликоль) в рамках единого подхода. 

При этом критериями поиска перспективных стабилизаторов для кластеров меди будем считать синтетиче- 
скую доступность, возможность модификации заместителей для создания максимальных стерических препятствий к 
поверхности кластера металла, свойства восстановителя и формирование наиболее прочных связей с медью (макси- 
мальная энергия комплексообразования). 

Экспериментальная часть. В рамках исследования были рассмотрены кластеры Си, (п=1-7, 13). Для иссле- 
дования влияния электронного строения азотного центра на процесс комплексообразования были использованы: про- 
изводные аммиака — карбамовая кислота Н›М-СООН и анилин Н>М-РВ; молекулы с азотным центром в системе со- 
пряжения — пиррол НМС.Н. и пиридин МС5Н5; структуры азота с двойной связью — НМ=СН) (метанимин), НМ=МН 
(диазен) и НМ=О (нитроксил); атомы азота с тройной связью — тиоцианат лития МС$“Та, цианид лития М=С”ТА и азот 
№. Для сравнения рассчитаны некоторые комплексы с молекулярным кислородом и водой. Расчеты проводились с 
полной оптимизацией всех параметров методом теории функционала плотности РВЕРВЕЛ ап! 20%, [13—15] с использо- 
ванием пакета программ Саиз1ап-09 [16]. Идентификация стационарных точек осуществлялась путем расчета матри- 
цы силовых постоянных. Все приведенные далее структуры являются минимумами на соответствующих поверхностях 
потенциальной энергии (ШПТЭ). Энергия стабилизации комплексов, образованных кластерами металла оценивалась 
как разность полной энергии комплекса и суммы полных энергий индивидуальных молекул. При наличии нескольких 
изомеров комплексов в таблицу вносился результат для наиболее стабильной формы. Влияние среды учитывалось в 
неявном виде — при помощи процедуры ЗСВЕ. Для исследования влияния среды выбраны апротонные среды различ- 
ной полярности — н-гептан и этиленгликоль, а также вакуум. 

Обсуждение и результаты. Расчеты методом РВЕРВЕЛ аШ207, показывают, что для кластеров Си, (п=3—6) 
наибольшей стабильностью отличаются плоские изомеры, для которых основным структурным мотивом являются 
треугольные формы. В тоже время кластеры где п>6 стабилизируются лишь в ЗО формах. Так, на рис. 1 приведены 
три изомера Сиз, соответствующие минимуму на ППЭ и имеющие неплоское строение. 

Как показывают расчеты, результаты которых приведены в табл. | и на рис. 2, в газовой фазе все производные 
аммиака (карбамовая кислота и анилин) проигрывают по устойчивости своих комплексов с поверхностью кластеров 
меди, как воде, так и молекуле кислорода (-21-22ккал/моль). При этом, с ростом размера кластера, энергия стабили- 
зации таких структур изменяется нелинейно, подобно синусоиде. 
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АЕ=0 АЕ=5,21 АЕ=22,33 
Рис. 1. Геометрические характеристики минимумов, полученных при помощи расчетов 
методом РВЕРВЕЛ ап1207, для структуры Си7 
На рис. 1 все длины приведены в ангстремах, ДЕ — энергия относительной устойчивости в ккал/моль. Нуле- 


вое значение соответствует самому устойчивому изомеру. 
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Таблица 1 
Энергии стабилизации комплексов, содержащих связь Си-М, рассчитанные методом 
РВЕРВЕЛ ап1201, (в ккал/моль) в условиях газовой фазы 
Молекула Си Си> Сиз Си Си5 Сиб Си Си1з 
Н›М-СООН 7,29 13,54 17,78 17,15 
о 13,76 20,87 16,11 24,80 
МН, 
() 12,65 11,96 16,26 15,11 
Н 
№ 
| р 21,37 30,00 36,71 33,74 26,84 
М 
НМ=СН, 24,40 32,08 40,24 36,55 28,66 27,13 27,89 
НМ=МН 31,02 33,16 46,05 40,18 31,75 27,28 30,57 36,71 
НМ=О 46,92 38,04 58,89 48,90 43,37 33,07 39,47 52,28 
МСТ 16,82 32,81 42,65 39,27 31,58 
М=М 14,25 21,01 27,53 25,23 16,62 
МС 23,36 33,83 47,06 37,15 31,04 
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Рис. 2. Геометрические характеристики минимумов, полученных при помощи расчетов методом РВЕРВЕ/Гап!207, 
для комплексов карбамовой кислоты, анилина, пиридина и пиррола с кластером Сиз 


Переход к молекулам пиррола и пиридина, в которых атом азота включен в цепь сопряжения, повышает 
устойчивость получающихся комплексов с Си-—М№-взаимодействием. Так, молекула пиридина связывается с поверхно- 
стью медного кластера прочнее, чем молекулы воды или кислорода, и вполне может препятствовать их контакту с 
металлом, учитывая же существенный размер молекулы пиридина (радиус онзагеровской сферы для нее составляет 


3,68А), ее локализация на одном атоме меди защищает еще несколько ближайших атомов кластера. 
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Наиболее стабильные комплексы на поверхности медного кластера дают молекулы, в которых азотный центр 
связан двойной связью с гетероатомом, при этом сила взаимодействия «кластер — молекула» растет с увеличением 
электроотрицательности заместителя в ряду НМ=СН, < НМ=МН < НМ=0. Некоторые из этих комплексов изображены 
на рис. 3. Следует отметить, что стабильность образованных структур значительно превышает силу контакта с по- 
верхностью меди молекул кислорода или воды (на -—10-20 ккал/моль). Как и в случае молекулы пиридина, энергия 
контакта Си-—М изменяется нелинейно, возрастая до максимума в случае кластера Сиз и далее уменьшаясь до миниму- 
ма в случае комплексов с кластером Си. Следующая зона роста энергии Си_—М№-взаимодействия начинается с семи- 
атомного кластера меди (п=7). Для икосаэра Си: энергия Си-М взаимодействия сопоставима по величине с таковой 
для трехатомных кластеров меди, при этом координация с поверхностью медного кластера идет исключительно по 


азотному центру молекулы, без участия непроделенных электронных пар заместителей. 
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Рис. 3. Геометрические характеристики минимумов, полученных при помощи расчетов 


методом РВЕРВЕЛап!207, для комплексов метанимина с кластерами Си, (п=1"6) 


Органические молекулы, в которых азот связан тройной связью, занимают промежуточное положение по силе 
стабилизации их комплексов с кластерами меди, при этом структуры Си, “М№5 с атмосферным азотом по всем пунктам 
проигрывает комплексам Си„”О., а токсичные соли цианида и тиоцианата лития показывают результаты близкие к 
ранее рассмотренным молекулам НМ№=У (У=СНЬ, МН, О), но более низкие по энергетике. Некоторые такого рода 
структуры приведены на рис. 4. Недостатком таких систем является их токсичность и возможность модификации 
только за счет противоиона, что дает немного вариантов, мало различающихся между собой. 

Таким образом молекулы, в которых азот связан двойной связью с кислородом в условиях газовой фазы, яв- 
ляются перспективными структурными мотивами сложных молекулярных стабилизаторов кластеров меди, при этом 
токсичность молекул, полученных на их основе, будет достаточно низкой. 
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Рис. 4. Геометрические характеристики минимумов, полученных при помощи расчетов 
методом РВЕРВЕ/Лап!207, для комплексов азота, цианида и тиоцианата лития с кластером Сиз 


Переход к неполярной среде п-гептана ослабляет все Си_—М-взаимодействия (табл. 2), но ослабление по- 
разному зависит как от размера кластера, так и от природы азотсодержащей молекулы. При этом все обнаруженные 
ранее тенденции сохраняются и наибольшей стабильностью отличаются структуры, где азотный центр связан двойной 
связью с гетероатомом, при этом максимальную эффективность по-прежнему показывает молекула нитроксила 
НМ=О. Хотя молекулы диазена НМ=МН и метанимина НМ=СН»› ощутимо проигрывают ей в качестве молекулы- 
стабилизатора в уловиях среды предельных углеводородов, они также являются перспективными моделями для поис- 
ков присадок за счет большего количества центров для замещения. В условиях неполярной среды соли цианида лития 
и тиоцианата лития заметно ухудшают свои свойства по сравнению с молекулами НМ=У (У=СНь, МН, О). 



































Таблица 2 
Энергии стабилизации комплексов, содержащих связь Си-М, рассчитанные методом 
РВЕРВЕ/Тап1207, (в ккал/моль) в неполярной системе (н-гептан) 
Молекула Си! Си> Сиз Си Сиз Сиб 
Н›М-СООН 7,33 13,70 17,63 16,13 
13,97 21,33 25,92 23,93 
МН, 
() 12,43 11,75 15,55 13,61 
Н 
№ 
| р. 19,64 29,88 35,54 32,35 
М 
НМ=СН» 17,42 33,27 37,79 35,41 29,16 24,45 
НМ=МН 28,82 30,39 43,31 36,48 32,04 24,57 
НМ=О 46,85 35,52 58,35 46,30 45,05 
МСТ 10,92 25,23 33,81 30,39 
№ 14,43 20,71 27,03 24,06 
МСТ 17,42 33,27 37,79 35,41 31,34 29,69 





























Изменение среды на полярную апротонную среду двуатомного спирта (этиленгликоль) приводит к усилению 
всех Си—М-взаимодействий на поверхности медных кластеров с сохранением всех обнаруженных ранее закономерно- 
стей и тенденций как по сравнению с газовой фазой, так и по сравнению со средой предельного углеводорода. Все 
соответствующие результаты приведены в табл. 3. 
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Как показывают расчеты, в полярной среде многоатомных спиртов производные аммиака (карбамовая кисло- 
та и анилин) существенно проигрывают в силе образованных ими Си—М№-комплексов азотсодержащим молекулами с 
эффектами сопряжения (пиридин, пиррол) и молекулам с тройными связями (тиоцианаты и цианиды щелочных ме- 
таллов). При этом наиболее прочные связи по-прежнему регистрируются для молекул с двойными связями, с без- 
условным лидером — молекулой нитроксила МН=О. 



































Таблица 3 
Энергии стабилизации комплексов, содержащих связь Си-М, рассчитанные методом 
РВЕРВЕ/Л ап1201, (в ккал/моль) в полярной апротонной системе (этиленгликоль) 
Молекула Си Си? Сиз Си Си5 Сиб 
о 14,61 21,79 25,66 24,44 
МН, 
® 12,47 11,37 14,49 13,38 
Н 
< 
| В 21,44 | 29,28 33,78 32,09 | 26,7 
М 
НМ=СН, 23,50 31,59 37,01 34,89 28,65 27,22 
НМ=МН 32,75 31,92 46,13 38,82 31,03 26,58 
НМ=О 54,07 36,70 65,00 50,25 47,62 33,11 
МСУ 12,82 26,76 25,67 
№ 14,88 19,91 24,26 
МЕС 32,36 31,13 34,93 29,65 28,53 























Выводы. Рассмотрев все основные случаи простейших азотсодержащих молекул, перспективных для исполь- 
зования как стабилизаторов нанопорошка меди, авторы выбрали для дальнейших экспериментов нитроксил Н-М=О, 
диазен Н-М=М-Н и метанимин Н-М=СН», которые являются синтетически доступными, прочно связывающимися с 
поверхностью кластера меди молекулами, также проявляющими восстановительные свойства и не нарушающими 
структуры кластера меди при координации. При этом, заменяя водородный атом на любой объемный заместитель, 
возможно блокировать доступ к большой площади поверхности кластера различных агрессивных соединений. Следу- 
ет также отметить, что комплексы кислорода на поверхности медного кластера существенно уступают в прочности 
структурам, образованным молекулами НМ=У (У=СНь», МН, О). Полученные данные позволяют оптимизировать син- 
тез стабилизированных медных наночастиц для их дальнейшего применения в качестве функциональных наноматери- 
алов триботехнического назначения. 
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